Stripe Phases in Layered Nickelates 
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To describe quasi two-dimensional nickelates we introduce an efltective Hamiltonian for Cg electrons 
which includes the kinetic energy, on-site Coulomb interactions, spin-spin and Jahn- Teller (static) 
terms. The experimental stripe phases are correctly reproduced by the model. The mechanisms 
responsible for stripe formation are different than those reported in cuprates and manganites. 
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Strong electron correlations are responsible for numer- 
ous interesting properties of doped transition metal ox- 
ides. In these systems coexisting charge, magnetic and 
orbital order emerge from competition of different type 
of electronic energy and the coupling with the lattice 
Among them the phenomenon of stripes is common 
— it plays an important role in cuprates j2|] and occurs 
also in nickelates [3] and manganites [1]. The stripes 
in cuprates arise from the competition between kinetic 
energy of doped holes and magnetic energy of ordered 
spins • The evolution of metallic stripes under increas- 
ing hole doping and their spectral properties could be 
described within a purely electronic model Q. 

In contrast, the origin of stripe phases which involve 
orbital order is more subtle. The kinetic energy is here 
anisotropic and the orbital states easily couple to the 
lattice. The extension by Jahn- Teller (JT) interactions 
is necessary in the models describing doped manganites 
and nickelates ['^J. The stripes in monolayer nicke- 
lates La2-xSri+2a:Ni04 were observed in a range of dop- 
ing 0.289 < x < 0.5 3]. They occur as charge walls in an 
antiferromagnetic (AF) phase and were described within 
the electronic model in a broad range of doping x < 0.4 
[§| , while the electronic structure calculations reproduced 
them at X = 1/3 and 1/2 Here we describe stripe 

phases within an effective model featuring only eg elec- 
trons at Ni sites renormalized by oxygen ions. At each 
site the local basis is given by two Cg orbitals symmetry. 



I.e. 



y 



^ and z = 3z^ — orbitals. 



The degenerate Hubbard Hamiltonian H for two Cg 
orbital states contains several terms, 

"H = iJkin + Hci + Hint + -ffjT- (1) 

The kinetic part iJkin is given by (dda) element t, 

r 
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and contains the hopping between x and z orbitals at 
neighboring sites {ij} which changes sign between bonds 



along a and h axis. The creation operators d\^^ corre- 
spond to electron in orbital fi = x,z, with spin a =t,i 
located at the site i. The kinetic energy is supplemented 
by orbital splitting > due to crystal field 



r) = -E, 



(3) 



iJint stands for on-site Coulomb interactions between eg 
electrons and was used before for monolayer, bilayer and 
cubic manganites [i!] and does not need to be reproduced 
here; it is parametrized by intraorbital Coulomb element 
U and Hund's exchange Jh- 

The simplified JT interaction iJjT is: 



H 



JT 
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where = Y^aa 
the pseudospin r = 



ta is eg electron density operator and 
1/2 operators are defined as follows: 



2 '^^XAxa'^izcy 
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The JT term includes three modes, {Qu, Q2i, Qsi} which 
denote standard eg static deformations of the octahedron 
around site i. Note that for dynamic JT effect each two 
neighboring Ni^+ ions share one oxygen and thus the JT 
distortions around them are not independent. 

We performed Hartree-Fock (HF) computations on a 
finite 6x6 cluster (with periodic boundary conditions) 
filled by N electrons and determined the ground state by 
comparing energies of several converged HF states. Fol- 
lowing them the electron correlations were studied using 
Local Ansatz (LA) method which implements the lead- 
ing on-site electron correlations [11| . Due to high density 



of eg electrons per site n = 2 — x the correlation energies 
were found to be much larger here than those reported 
before for cuprates 12 1 or manganites @, but they do 



not modify the stability of the HF ground state. 
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FIG. 1: The antiferromagnetic ground state obtained for 
La2Ni04 (at a; = 0). Each circle represents Cg electron den- 
sity with its diameter corresponding to one-half of density; 
the arrow length — one-half of the eg spin; the horizontal 
bar length (here absent) — charge density difference between 
X and z orbitals (longest bar to the right/left corresponds 
to pure x/z orbital state). Parameters (in eV): t = 0.6, 
Jh = 0.9, E, = -0.6; U = 8to, p.jT = 3.0 eV/A, K = eV/A^ 



In the undoped La2Ni04 (TV = 2x6^ corresponds 
to n = 2 electrons per site) we found a uniform charge 
distribution in the ground state, with equal orbital elec- 
tron densities n^; = = 1, corresponding to a high-spin 
S = 1 state at each Ni^+ ion. The ionic spins S* = 1 are 
coupled by the AF interaction which follows from the 
Hamiltonian % Eq. (1) in a similar way as it does for 
s = 1/2 spins in cuprates. This explains the origin of 
G-type AF order (G-AF) shown in Fig. 1. 

Next we investigate the ground states for decreasing 
electron number per site n — 2 — x. At low doping 
0.1 < X < 0.3 we obtained within the present approach 
isolated polaronic states (not shown). Diagonal stripe 
phases are found for doping x = l/i and x = 1/2, in 
agreement with experimental observations Q — they are 
shown in Figs. 2 and 3. They are not like the stripes in 
cuprates or manganites. In cuprates charge order coex- 
ists with the modulation of AF order between alternating 
AF domains, separated by nonmagnetic and half-filled 
domain walls. The JT coupling does not contribute here 
as the hole sites are inactive. In contrast to cuprates, 
Ni ions always carry magnetic moments and contribute 
to the magnetic order within a single domain of the G- 
AF phase. Note also that the present stripe phases are 
insulating, while the ones in cuprates are metallic [5| . 

In manganites the situation is more complex due to 
possible presence of orbital order and due to extra Hund's 
exchange coupling of Cg electrons with t2g core spins 
S = 3/2 1]. Here several competing mechanisms en- 
ter on equal footing and distinct stripe patterns with the 
true long range order are not a rule. As a result, the 
magnetic order of the ground state was found to be very 



OQ0QO© 
0©00®Q 
©QO0OQ 
OQ0Q©® 
0©00®Q 
©OO©©© 



©©© 
©®© 
©©© 
©©© 
©©© 
©©© 



©©© 
©©© 
©©© 
©©© 
©©© 
©©© 



FIG. 2; Two stripe phases with C-AF spin order (FM ver- 
tical and AF horizontal order) obtained for doped nickelates 
at doping a; = 1/3 (top) and a; = 1/3 (bottom). Upper panel 
— charge minority sites form diagonal (11) stripe bound- 
aries (each third line). Lower panel — Perfect checkerboard 
crystal-like order given by charge majority/minority sublat- 
tices. In both cases holes are doped predominantly into x 
orbitals. Legend and parameters as in Fig. 1. 



sensitive to the precise values of the model parameters 
13j . However, at x — 1/2 doping one finds a robust 
instability for monolayer manganites toward two sublat- 
tices with charge majority and minority sites, see also 
[8a] . The JT effect stabilizes here the CE magnetic order 
accompanied by checkerboard charge order, with the JT 
distortions contributing only on charge majority sites. 

In monolayer nickelates one finds the opposite situa- 
tion: The main difference in the ionic picture is that in 
the undoped substance there are two Cg electrons per site, 
and hole doping generates JT active sites, where a single 
electron couples to JT distortions Q2i and Qsij see Fig. 
3. This coupling is crucial for stabilizing ordered stripe 
structures at x = 1/3 and a; = 1/2 doping, shown in Fig. 
2. In effect one can say that the kinetic energy of doped 
holes is suppressed and the competition between mag- 
netic and kinetic energy along the charge minority walls 
is quenched, in agreement with double exchange [j]. 

The above ionic picture used to interrpret both stripe 
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sites can order at doping x = 1/3 and appear as diagonal 
stripes (each third hne). Similarly, they give diagonal 
stripes (each second line) for x = 1/2 doping stabilized 
by the JT distortions. Indeed, when the JT coupling is 
turned off (at (?jt = 0), the stripes do not form. Strong 
electron correlations are also important and we have ver- 
ified that similar stripe phases to those shown in Figs. 
2 and 3 occur for a stronger on-site Coulomb repulsion 
U = 12t in the present model Eq. (1), only the charge 
modulation is somewhat enhanced. 

In summary, we identified the coupling of Cg electrons 
at sites doped by holes to local Jahn- Teller distortions as 
the microscopic origin of diagonal stripe phases in mono- 
layer nickelates. The presented results elucidate funda- 
mental difference between stripe phases in doped nicke- 
lates, cupratcs and manganitcs, as the ions active for the 
Jahn- Teller effect play a different role in each case. 
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FIG. 3: Jahn- Teller distortions corresponding to the stripe 
phases shown in Fig. 2: top — x = 1/3, and bottom x = 1/2 
doping. The distortions Qa (in A, scaled by a factor of 2) axe 
shown as vertical bars drawn slightly to the left of each site 
and Qai, by bars to the right of each site (isotropic mode Qn 
is absent). Legend and the values of parameters as in Fig. 1. 

phases is not far from the actual charge distribution 
obtained in the present computations. We have found 
rii = 1.82, 712 = 1-37 (ni = 1.66, n2 = 1.34) at charge 
majority/minority sites at a; = 1/3 {x = 1/2) doping. As 
quantum fluctuations are suppressed in the HF method, 
this corresponds to magnetic moments mi = ±0.89, 
1712 = ±0.67 (mi = ±0.81, m2 = ±0.62) in both phases. 

Doping occurs predominantly in x orbitals, where also 
the magnetic moments are suppressed. The density in 
X orbitals varies between Ux = 0.47(0.85) and rix = 
0.45(0.73) for the stripe phases at x = 1/3 and a; = 1/2 
doping, while density in z orbitals is > 0.89 in all 
cases. For this charge distribution both Q2i and Qsi JT 
modes are active at charge minority sites, sec Fig. 3. 
We suggest that JT distortions govern the stripe forma- 
tion in monolayer nickelates, and stabilize the obtained 
charge alternation accompanied by large JT distortions 
at charge minority sites. We have found that the doped 
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